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ABSTRACT 

The near-infrared colors of the planets directly imaged around the A star 
HR 8799 are much redder than most field brown dwarfs of the same effective 
temperature. Previous theoretical studies of these objects have concluded that 
the atmospheres of planets b, c, and d are unusually cloudy or have unusual 
cloud properties. Some studies have also found that the inferred radii of some or 
all of the planets disagree with expectations of standard giant planet evolution 
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models. Here we compare the available data to the predictions of our own set 
of atmospheric and evolution models that have been extensively tested against 
observations of field L and T dwarfs, including the reddest L dwarfs. We require 
mutually consistent choices for effective temperature, gravity, cloud properties, 
and planetary radius. This procedure thus yields plausible values for the masses, 
effective temperatures, and cloud properties of all three planets. We find that 
the cloud properties of the HR 8799 planets are not unusual but rather follow 
previously recognized trends, including a gravity dependence on the temperature 
of the L to T spectral transition-some reasons for which we discuss. We find 
the inferred mass of planet b is highly sensitive to whether or not we include 
the H and K band spectrum in our analysis. Solutions for planets c and d are 
consistent with the generally accepted constraints on the age of the primary star 
and orbital dynamics. We also confirm that, like in L and T dwarfs and solar 
system giant planets, non-equilibrium chemistry driven by atmospheric mixing is 
also important for these objects. Given the preponderance of data suggesting that 
the L to T spectral type transition is gravity dependent, we present an exploratory 
evolution calculation that accounts for this effect. Finally we recompute the the 
bolometric luminosity of all three planets. 

Subject headings: brown dwarfs — planetary systems — stars: atmospheres - 
stars: low mass, brown dwarfs - stars: individual (HR 8799) 



1. INTRODUCTION 



Establishing the masses, radii, effective temperatures, and atmospheric composition 
of the planets orbiting the A star HR 8799 has been a challenge. Of the four planets 



(Marois et al. 2008, 2010) directly imaged orbiting the star HR 8799, broad photometric 
coverage (1 - 5 jum) is available for three planets, b, c, and d (Marois et al. 2008 Currie 



et al. 2011), and some spectral data is available for one planet, b (Barman et al. 2011a). 



Efforts to fit the available data with atmosphere and evolution models have produced mixed 
results. In some cases the best-fitting models predict radii and ages that are at odds with 
other constraints, such as evolution models and the age of the system. The purportedly 
unusual cloud properties of the planets have also received great attention. Here we present 
an examination of the properties of HR 8799 b, c, and d using publicly available data as 
well as our own evolution and atmosphere models. Our aim is to determine if a set of 
planet properties can be derived that simultaneously satisfy all observational and theoretical 
constraints and to ascertain the nature of atmospheric condensate layers in each planet. 
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We open below with a summary of the model parameters previously derived for these 
planets. In the remainder of this section we briefly review what is known about the atmo- 
spheric evolution of brown dwarfs and discuss the issues that have arisen to date in the study 
of the HR 8799 planets, particularly regarding the inferred cloud properties and planet radii. 
In succeeding sections we explore the nature of clouds in low-mass objects more deeply and 
present model solutions for the masses, effective temperatures (T eff ), and cloud properties of 
the planets. We find, as have all other previous studies, that clouds are present in the vis- 
ible atmosphere of these planets at lower effective temperatures than in typical field brown 
dwarfs. In agreement with Barman et al. (2011a) but unlike most other previous studies 



e.g., Bowler et al. 2010 Currie et al. 2011 Madhusudhan et al. 2011) we find that the 



clouds of the HR 8799 planets are similar to those found in field L dwarfs. 



1.1. Masses and Radii of HR 8799 b,c, and d 



In the HR 8799 b, c, and d discovery paper, Marois et al. (2008) derived the mass and 



effective temperature of each object in two ways. In the first method they computed the 
luminosity of each object and compared that to theoretical cooling tracks for young giant 
planets given the constraint of their estimated age of the primary star. In the second method 



they fit atmosphere models derived using the PHOENIX code (Hauschildt et al. 1999) to 



the available six-band near-infrared photometry (1 to 2.5 /im) to constrain T c g and logg, the 
two most important tunable parameters of atmosphere models. Radii of each planet were 
derived by comparing the model emergent spectra with the observed photometry and known 
distance to the target. Notably only models that included the effects of refractory silicate 
and iron clouds were consistent with the data. However the radii estimated by this method 
were far smaller than expected for solar metallicity gas giant planets at such young ages. 

A number of followup studies presented new data new data and models in an attempt 



to better understand the planets. Barman et al. (2011a) fit a suite of models to the available 



photometry (but not the M band (Galicher et al. 2011) data) and H and K band spectra 
that they obtained for planet b. By comparing the integrated flux from their best fitting 
model atmosphere to the estimated bolometric luminosity of the planet, they found a small 
radius for the planet R ~ 0.75 Rj. Galicher et al. (2011) also fit the Barman atmosphere 
models to the photometry, including new M band data. They found somewhat higher 



gravity solutions than Barman et al. (2011a) but also required a small radius for planet 
b, approximately 70%-or about one-third the volume-expected from planetary evolution 
models. Such a large discrepancy is difficult to reconcile with our understanding of both 
giant planet evolution and the high pressure equation of state of hydrogen. Instead the most 
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straightforward interpretation is that the atmosphere models are not representative of the 
actual planetary atmosphere and Barman et al. suggest that higher metallicity models might 
provide a better fit and give more plausible radii. 



Likewise Bowler et al. (2010) selected the model spectra (from among the models of 



Hubeny fc Burrows| Q2007fl ; |Burrows et al] Q2006fl ; |Allard et alT] ( |2001p ) which best fit the 
available photometry for HR 8799b. Their best fitting spectra were quite warm, with T e g 
from 1300 to 1700 K and thus they required even smaller radii (~ 0.4 Rj) in order to meet 
the total luminosity constraint given the photometry available at that time. 



In contrast Currie et al. (2011) searched for the best fitting models while requiring that 



the planet radii either matched those predicted by a set of evolution models (Burrows et 



al. 1997) or were allowed to vary. They found that what they termed to be "standard" 



brown dwarf cloud models required unphysically small planet radii to fit the data. However 
their "thick cloud" models could fit the data shortward of 3 /im by employing planetary radii 
that were within about 10% of the usual evolution model prediction. As we note below, 
however, the "standard" cloud model has itself not been demonstrated to fit cloudy, late 
L-type dwarfs; thus this exercise does not necessarily imply the planets' clouds are "non- 
standard". Nevertheless they were able to fit much of the photometry with planetary radii 
consistent with evolution model predictions. 



Finally Madhusudhan et al. (2011) explored a set of models similar to those studied by 



Currie et al. with yet another cloud model but without the radius constraint. Their best fits 
are very similar to those of Currie et al. but with somewhat lower T e g. 

The characteristics of the planets as derived in the 2011 publications are summarized in 
Table 1. Not all authors report every parameter so some radii and ages are left blank. Note 
the diverse set of masses, radii, and effective temperatures derived by the various studies. 
Despite the variety some trends are clear: planet b consistently is found to have the lowest 
mass and effective temperature and its derived radius is almost always at odds with the 
expectation of evolution and interior models. 

We note that at very young ages the model radii of giant planets depends on the initial 



conditions of the evolutionary calculation (Stevenson 1982 Baraffe et al. 2002; Marley et 



al. 2007a; Spiegel & Burrows 2012). However at ages younger than several hundred million 



years the planetary radius is expected to be no smaller than about 1.1 times that of Jupiter 



regardless of the formation mechanism. Hence radii derived by Barman et al. (2011a) and 



Galicher et al. (2011) are not consistent with evolutionary calculations, regardless of the 



initial boundary conditions. Indeed the equation of state for gas giant planets, even ones 
enriched in heavy elements, preclude such radii. 
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1.2. Clouds 



1.2.1. Brown Dwarfs 

As a brown dwarf ages it radiates and cools. When it is warm, refractory condensates, 
including iron and various silicates, form clouds in the visible atmosphere. Over time the 
clouds become progressively thicker and more opaque, leading to ever redder near-infrared 
colors. As the dwarf cools the cloud decks are found at higher pressures, deeper in the atmo- 
sphere. Eventually the clouds disappear from the photosphere. Indeed the first two brown 



dwarfs to be discovered, GD 165B (Becklin & Zuckerman 1988) and Gl 229B (Nakajima et 



al. 1995), were ultimately understood to represent these two different end cases: the cloudy 
L and the clear T dwarfs (see Kirkpatrick (2005) for a review). Understanding the behavior 



of clouds in substellar atmospheres and how it might vary with gravity has become one of 
the central thrusts of brown dwarf science. 

The earliest models for these objects assumed that the condensates were uniformly 
distributed vertically throughout the atmosphere (e.g., Chabrier et al. 2000). Later, more 



sophisticated approaches attempted to model the formation of discrete cloud layers that 
would result from the gravitational settling of grains. 

With falling effective temperature, T e g, the bases of the iron and silicate cloud decks are 
found progressively deeper in the atmosphere. Because of grain settling the overlying atmo- 
sphere well above the cloud deck loses grain opacity and becomes progressively cooler. Thus 
over time more of the visible atmosphere becomes grain free and cooler. Cooler temperatures 
favor CH4 over CO. The removal of the opacity floor that the clouds provided at higher T e g 
also allows flux in the water window regions to escape from deeper in the atmosphere. This 
leads to a brightening in the J band and a blueward color shift in the near-infrared. In field 
brown dwarfs this color change begins around effective temperature T e g- ~ 1200 to 1400 K 
and is complete over a strikingly small effective temperature range of only 100 to 200 K (see 



Kirkpatrick (2005) for a review). This experience led to the presumption that all objects 
with effective temperatures below about 1100 K would have blue near- infrared colors, like 
the field brown dwarfs. 



1.2.2. HR 8799 b, c, and d 



The early directly imaged low mass companions confounded these expectations from the 
brown dwarf experience. The companion 2MASSW J1207334-393254 b (hereafter 2M1207 
b) has red infrared colors despite its low luminosity and apparently cool T e g ( Chauvin et al. 
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2004) . Likewise the HR 8799 planets have colors reminiscent of hot, cloudy L dwarfs but 



their bolometric luminosities coupled with radii from planetary structure calculations imply 
T cS ~ 1000 K or lower (iMarois et al. 1120081 |20Tob. 



The red colors, particularly of the HR 8799 planets, spawned a storm of studies inves- 
tigating the atmospheric structure of the planets. Essentially all of these papers concluded 
that the planets could be best explained by invoking thick cloud decks. Since this ran 
counter to expectation, these clouds were deemed "radically enhanced" when compared to 
"standard" models ( Bowler et al.|2010~ ). Likewise Currie et al. ( 2011[ ) compared their data to 



the Burrows et al. (2006) model sequence and concluded (their §5) that the HR 8799 planets 



have much thicker clouds than "...standard L/T dwarf atmosphere models." Madhusudhan 



et al. (2011) state that their fiducial models "...have been shown to provide good fits to 



observations of L and T dwarfs (Burrows et al. 2006)". They then find that much cloudier 



models are required to fit the imaged exoplanets and thus conclude that the cloud properties 
must be highly discrepant from those of the field L dwarfs. 

Such conclusions, however, seem to overlook that the study of L dwarf atmospheres 
is still in its youth. Cloudy atmospheres of all kinds are challenging to model and the L 
dwarfs have proven to be no exception. Thus whether or not the HR 8799 planets have 
unusual clouds depends on the point of reference. Indeed while most published models of 
brown dwarfs are able to reproduce the spectra of cloudy, early L-type dwarfs and cloudless 
T dwarfs, the latest, reddest — and presumably cloudiest — L dwarfs have been a challenge. 



The points of comparison for the work of Currie et al. (2011) and Madhusudhan et al. (2011) 



were the models described in Burrows et al. (2006). When compared to the red-optical and 



near-infrared photometry of L and T dwarfs, those models did not reproduce the colors of the 



latest L dwarfs as the models are too blue (see figure 17 of Burrows et al. (2006)) implying 



that they lacked sufficient clouds. Burrows et al. (2006) also presented comparisons of their 



models to L dwarf spectra; however the comparisons are only to an LI and an L5 dwarf. 
There are no comparisons to very cloudy late L dwarf spectra in the paper so the fidelity 
of their model under such conditions cannot be judged. For these reasons a comparison 
of the cloudy HR 8799 planets to the "standard" L dwarf models, such as presented by 



Madhusudhan et al. (2011) and Currie et al. (2011), does not address the question whether 



the HR 8799 planets are really all that different from the cloudiest late L dwarfs since those 
models have apparently do not reproduce the colors of the latest L dwarfs. 

At least one set of atmosphere and evolution models is available that has been compared 



against the near- to mid-infrared spectra and colors of latest L dwarfs. In Cushing et al. 



(2008) and Stephens et al. (2009) we compared our group's models to observed far-red to 



mid-infrared spectra of L and T dwarfs, including L dwarfs with IR spectral types as late as 
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L9 (with 7 objects in the range L7 to L9.5). We found that the models with our usual cloud 
prescription fit the spectra of L dwarfs of all spectral classes (including the latest field dwarfs) 
well, but not perfectly. In Saumon & Marley (2008) we also presented a model of brown 
dwarf evolution that well reproduced the usual near-infrared color magnitude diagrams of 
L and T dwarfs, including the reddest L dwarfs. Here we apply our set of cloudy evolution 
models to the HR 8799 planet observations in an attempt to better understand these objects. 



1.3. Chemical Mixing 



Shortly after the discovery of Gl 229B, Fegley & Lodders (1996) predicted that — as in 
Jupiter — vertical mixing might cause CO to be overabundant compared to CH4 in chemical 
equilibrium in this object. This was promptly confirmed by the detection of CO absorption 



at 4.6 /iin by Noll et al. (1997) and Oppenheimer et al. (1998). The overabundance is caused 



by the slow conversions of CO to CH 4 relative to the mixing time scale. 

An obvious mechanism for vertical mixing in an atmosphere is convection. Brown dwarf 
atmospheres are convective at depth where the mixing time scale is short (minutes). The 
overlying radiative zone is usually considered quiescent but a variety of processes can cause 
vertical mixing, albeit on much longer time scales. Since the conversion time scales for 
CO — > CH4 and N2 — > NH3 range from seconds (at T ~ 3000 K) to many Hubble times (for 
T < 1000 K), even very slow mixing in the radiative zone can drive the chemistry of carbon 
and nitrogen out of equilibrium. From this basic consideration, it appears that departures 
from equilibrium are inevitable in the atmospheres of cool brown dwarfs and indeed the 



phenomenon is well established (e.g., Saumon et al. 2000; Geballe et al. 2001 Hubeny & 



Burrows|[2007l |Geballe et a"Ll[2009l |Mainzer et aLl[2007l |Saumon et afl|2006l |Stephens et~aL 



2009). 



With falling gravity the point at which chemical reactions are quenched occurs deeper in 
the atmosphere, where the higher temperature result in a greater atmospheric abundance of 
CO ( Hubeny fc Burrows|2007 Barman et al.||2011a ). At exoplanet gravities, mixing can even 



produce CO/CH4 ratios in excess of 1 (Barman et al. 2011a). Thus a complete giant planet 
exoplanet atmosphere model must account for such departures from chemical equilibrium as 
well. 
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2. Gravity, Refractory Clouds and the L/T Transition 

2.1. Nature of the Transition 

Two main causes of the loss of cloud opacity at the L to T transition have been suggested. 
In one view the atmospheric dynamical state changes, resulting in larger particle sizes that 
fall out of the atmosphere more rapidly, leading to a sudden clearing or collapse of the cloud 
(Knapp et al. 2004 Tsuji &; Nakajima 2003; Tsuji et al. 2004). This view is supported by 



fits of spectra to model spectra (Saumon & Marley 2008) computed with the Ackerman & 



Marley (2001) cloud model. In that formalism, a tunable parameter, / se <j controls cloud 
particle sizes and optical depth. Larger / se d yields larger particles along with physically and 
optically thinner clouds. Cushing et al.| (2008) and Stephens et al. (2009) have demonstrated 
that progressively later dwarfs (L9 to T4) can be fit by increasing / sed across the transition 
at a nearly fixed effective temperature. A variation on this hypothesis is that a cloud particle 



size change is responsible for the transition (Burrows et al. 2006). 



The second view is inspired by thermal infrared images of the atmospheres of Jupiter 
and Saturn at ~ 5um U 



Westphal 


1969 


Westphal et al. 


1974 


Orton et al. 


1996 


Baines 



et al. 2005). Gaseous opacity is low at this wavelength and the clouds stand out as dark, 
mottled features against a bright background of flux emitted from deeper, warmer levels 
in the atmosphere. Such images of both Jupiter and Saturn clearly show that the global 



cloud decks are not homogenous, but rather are quite patchy. Ackerman & Marley (2001) 



Burgasser et al. (2002), and Marley et al. (2010) have suggested that the arrival of holes 



in brown dwarf clouds, perhaps due to the clouds passing through a dynamical boundary in 
the atmosphere, might also be responsible for the L to T transition. This view is supported 
by the discovery of L-T transition dwarfs that vary in brightness with time with relatively 
large near- infrared amplitudes ( Artigau et al.|2009~ Radigan et al.||2011 ). Indeed Radigan (in 
prep) has found in a survey of about 60 L and T type brown dwarfs that the most variable 
dwarfs are the early T's, which are in the midst of the J — K color change. 

In order to match observations, modern thermal evolution models for the cooling of 
brown dwarfs have to impose some arbitrary mechanism, such as varying sedimentation 
efficiency or the imposition of cloud holes, by which the thick clouds in the late L dwarfs 
dissipate. A uniform cloud layer that simply sinks with falling T e g as the atmosphere cools 
turns to the blue much more slowly than is observed. Application of such a transition 
mechanism to reliably reproduce the colors and spectra of late L and early T dwarfs (e.g., 
near-infrared color-magnitude diagrams) led to the expectation that the normal behavior for 
cooling brown dwarfs-or extrasolar giant planets-is to turn to the blue at around 1300 K. 



However there have been indications that such a narrative is too simplistic and that 
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gravity plays a role as well. Two brown dwarf companions to young main sequence stars were 
found to have unexpectedly cool effective temperatures for their L-T transition spectral types 
by Metchev & Hillenbrand (2006) and Luhman et al. (2007). The analysis of Luhman et al. 



of the T dwarf HN Peg B was further supported by additional modeling presented in Leggett 



et al. (2008). Dupuy et al. (2009) presented evidence of a gravity dependent transition T e g 



on the basis of a dynamical mass determination of an M8 + L7 binary. Stephens et al. (2009) 
fit the model spectra of Marley et al. (2002) to the 1-15 /zm spectra of L and T dwarfs and 
found that L dwarf cloud clearing (as characterized by large f se d) occurs at T c g ~ 1300 K for 
log g = 5.0 and at ~ HOOK for log g = 4.5, although the sample size was admittedly small 
(Figure 1). Nevertheless such an association implies a cooler transition temperature at even 
lower gravity. 



2.2. Clouds at Low Gravity 

Even if directly imaged planets are not considered, there is already considerable evi- 
dence that the cloud clearing associated with the L to T transition occurs at lower effective 
temperatures in lower gravity objects than in high gravity ones. To understand what un- 
derlies this trend it is necessary to consider three separate questions. First, where does 
the optically-thick portion of the cloud lie in the atmosphere relative to the photosphere, 
as a function of gravity? An optically-thick cloud lying well below the photosphere will be 
essentially invisible whereas the same cloud lying higher in the atmosphere would be easily 
detected. Second, how does the total optical depth of the cloud vary with gravity? This is a 
complex problem involving the pressure of the cloud base and the particle size distribution. 
Third, how does the mechanism by which clouds dissipate vary with gravity? For example, 
do holes form at a different effective temperature in different gravity objects? In this section 
we consider only the first two questions and defer the third question to Section 5.6. 

To address the first question we need to understand how atmospheric temperature T 
varies with pressure P as a function of gravity. For a fixed effective temperature, a lower 
gravity atmosphere is warmer at a fixed pressure level than a higher gravity one. This is 
because more atmospheric mass-and thus greater opacity-overlies a given pressure level at 
lower gravity. Figure 2 provides an example using our model profiles. Since at equilibrium 
condensation begins at the intersection of the vapor pressure and thermal profiles, the cloud 
base occurs at lower pressure (higher in the atmosphere) in a low gravity object than a high 
gravity one. 

As objects cool with time (at essentially fixed gravity) clouds will persist at lower pres- 
sure and remain visible to cooler effective temperatures in lower gravity objects than higher 
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gravity ones. For example in Figure 2 the lowest gravity model shown at T e g = 900 K is hot- 
ter at all pressures greater than a few hundred millibar than a higher gravity T c g = 1300 K 
object. As explained below this degeneracy between cooler low gravity and warmer high 
gravity temperature profiles lies at the heart of the problem of simultaneously distinguishing 
gravity and effective temperature with a limited photometric dataset. 

Addressing the second question requires us to understand how the cloud column optical 
depth varies with gravity. This depends both on the amount of condensible material in the 
atmosphere available to form clouds and on the cloud particle size. From basic scaling laws 



and mass balance Marley (2000) derived an expression for the wavelength-dependent total 



column optical depth r A of a cloud in a hydrostatic atmosphere 

n = l 10 ^^) C-^) C-^^- ) . (1) 



lbar/ V g ) \ r eff 

Here P c i, r eff and p c refer to the pressure at the cloud base and the condensate effective 
(area-weighted) radius^ and density (see also Eq. 18 of Ackerman & Marley (2001)). <p 



is the product of the condensing species number mixing ratio and the ratio of the mean 
molecular weight of the condensate to that of the atmosphere. The expression assumes that 
some fraction e of the available mass above the cloud base forms particles with extinction 



cross section Q\ (which can be computed through Mie theory) . Ackerman & Marley (2001) 
also estimate the column optical depth of a cloud with a similar result. Generalizing their 
Eq. 16, 

gr cS (l + / sed ) 

Both Equations (1) and (2) hold that all else being equal-including particle sizes-we 
expect r oc P c i/g, just because the column mass above a fixed pressure level is greater at 
low gravity and there is more material to condense. Any cloud model which self-consistently 
computes the column mass of condensed material should reproduce this result. As shown 
above, however, the cloud base is at lower pressure in lower gravity objects, roughly P c i oc g, 
thus leading to the expectation that the cloud r would be approximately constant with 
changing gravity. This is not exactly true since there is a slope to the vapor pressure 
equilibrium curve and thus the actual variation is somewhat weaker, but the effects of gravity 
and the cloud base pressure alone do not strongly influence cloud column optical depth. 

The second component affecting the column cloud opacity is particle size. While a cloud 
model is required for rigorous particle size computation, we can examine the scaling of size 



1 Marley (2000) employed the mean particle size r c rather than the more rigorous area-weighted size. 
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with gravity. At lower gravity particle fall speeds are reduced, which reduces the downward 
mass flux carried by condensates of a given size r. Since fall speed is proportional to r 2 in the 
Stokes limit (the viscous regime at low Reynolds numbers) while the mass is proportional 
to r 3 , the flux scales with r 5 , a slight increase in particle size can produce the same mass 
balance in the atmosphere at lower gravity, and thus r is expected to increase relatively 
slowly with decreasing g. At large Reynolds number the dependence on fall speed is weaker 



than r and the equivalent result is found. Indeed recasting the Ackerman & Marley (2001) 
model equations suggests r oc (fsed/g) 1 ^ 2 , although the actual dependence is more complex 
as it depends upon an integral over the size distribution. Tests with the complete cloud 
model coupled to our atmosphere code predict about a factor of 4 increase in cloud particle 
radius (25 to 100 /xm) as gravity decreases by an order of magnitude from 300 to 30 ms -2 , a 
slightly faster increase than Jg. A roughly r oc g^ 1 ^ 2 relationship is also seen in the cloud 



model of Cooper et al. (2003) (see their Figures 2, 3, and 4). Returning to Eq. (1) and 
combining with the scaling discussed above thus suggests that all else being equal we expect 
cloud t oc Jg. 

Figure 3 illustrates all of these effects in model cloud profiles calculated for three atmo- 
sphere models with varying g and T e g. The atmospheric gravity spans two orders of magni- 
tude while the effective temperature varies from 1200 to 1000 K from the warmest to coolest 
object. As expected the cloud particle size indeed varies inversely with gravity(r ~ g^ 1 ^ 2 ) 
while the cloud base pressure decreases with decreasing gravity. The choice in the plot of a 
cooler T e ff for the lowest gravity object counteracts what would otherwise be an even greater 
difference in the cloud base pressure. The net result is that the total column optical depth 
for the silicate cloud in all three objects is very similar, r ~ 10. Thus a cooler, low gravity 
object has a cloud with a column optical depth that is almost indistinguishable from that of 
a warmer, more massive object. 

The thicker portion of the lines denoting cloud column optical depth signify the regions 
in the atmosphere where the brightness temperatures between A = 1 and 6 /im are equal to 
the local temperature. In other words the thick line represents the near-infrared photosphere. 
In all three cases there is substantial cloud optical depth (r A > 0.1) in the deeper atmospheric 
regions from which flux emerges in the near-infrared. As a result clouds play comparable 
roles in all three objects despite the two order of magnitude difference in gravity and the 
200 K temperature difference. We thus conclude that the net effect of all of these terms is to 
produce clouds in lower gravity objects with optical depths and physical locations relative 
to the photosphere comparable to clouds in objects with higher gravity and higher effective 
temperature. 
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3. Modeling Approach 



To model the atmospheres and evolution of exoplanets we apply our usual modeling 
approach which we briefly summarize in this section. We stress that the fidelity of model fits 
in previous applications of our method to both cloudy and clear atmosphere brown dwarfs 



(Marley et al.|1996, 


2002 


Burrows et al.||1997 


Roellig et al.||2004 


Saumon et al.|2006, 


2007 


Leggett et al. 2007a|b; Mainzer et al. 2007; Blake et al. 2007 Cushing et al. 2008 Geballe 


et al. 


2009 




Stephens et al. 


2009 


) validates our overall approach and provides a basis of 



comparison to the directly imaged planet analysis. In addition to brown dwarfs the model 
has been applied to Uranus (Marley & McKay 1999) and Titan (McKay et al.|[l989) as well. 



3.1. Atmosphere and Cloud Models 



The atmospheric structure calculation is described in McKay et al. (1989); Marley et 



al. (1996); Burrows et al. (1997); Marley & McKay (1999); Marley et al. (2002); Saumon & 



Marley (2008). Briefly we solve for a radiative-convective equilibrium thermal profile that 



carries thermal flux given by crT e 4 ff given a specified gravity and atmospheric composition. 



The thermal radiative transfer follows the source function technique of Toon et al. (1989) 



allowing inclusion of arbitrary Mie scattering particles in the opacity of each layer. Our 



opacity database includes all important absorbers and is described in Freedman et al. (2008). 



There are, however, two particularly important updates to our opacity database since 



Freedman et al. (2008). First we use a new molecular line list for ammonia (Yurchenko et al. 



2011). Secondly we have updated our previous treatment of pressure-induced opacity arising 



from collisions of H 2 molecules with H 2 and He. This new opacity is discussed in Frommhold 



et al. (2010) and the impact on our model spectra and photometry in general is discussed in 



Saumon et al. (2012) 



The abundances of molecular, atomic, and ionic species are computed for chemical equi- 
librium as a function of temperature, pressure, and metallicity following |Fegley fc Lodders 



(1994 


1996 


); 


Lodders 


(1999 


); 


Lodders & Fegley 


(2002 


); 


Lodders 


(2003 


); 


Lodders & Fegley 


(2006 


assuming the elemental abundances of 


Lodders 


(2003 


). In this paper we explore only 



solar composition models. 



(2001) which 



For cloud modeling we employ the approach of Ackerman & Marley 
parameterizes the importance of sedimentation relative to upwards mixing of cloud particles 
through an efficiency factor, / se d. Large values of / se d correspond to rapid particle growth 
and large mean particle sizes. Under such conditions condensates quickly fall out of the 
atmosphere, leading to physically and optically thinner clouds. In the case of small / se d 
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particles grow more slowly resulting in a larger atmospheric condensate load and thicker 
clouds. Both our cloud model and chemical equilibrium calculations are fully coupled with 
the radiative transfer and the (P, T) structure of the model during the calculation of a model 
so that they are fully consistent when convergence is obtained. 

We note in passing that the cloud models employed in previous studies of the HR 8799 
planets have been ad hoc, as straightforwardly discussed in those papers. Particle sizes, 
cloud heights, and other cloud properties are fixed at given values while gravity, T e g, and 
other model parameters are varied. The methodology used here is distinct since in each 
case we compute a consistent set of cloud properties given a specific modeling approach, the 
Ackerman & Marley cloud. 

The coupled cloud and atmosphere models have been widely compared to spectra and 
photometry of L and T dwarfs in the publications cited in the introduction to this section. 



We emphasize in particular that Cushing et al. (2008) and Stephens et al. (2009) show 



generally good fits between our model spectra and observations of cloudy L dwarfs. The 
near-infrared colors of brown dwarfs are quite sensitive to the choice of / se( j, a point we will 
return to in Section 5.4. 



3.2. Evolution Model 



Our evolution model is described in Saumon & Marley (2008). In fitting the HR 8799 



data, we use the sequence computed with a surface boundary condition extracted from our 
cloudy model atmospheres with / se( j = 2. As we will see below, our best fits show that all 
three planets are cloudy with / se d = 2, which justifies this choice of evolution a posteriori. As 
the three planets appear to have significant cloud decks (as will be confirmed below), it is not 
necessary to use evolution sequences that take into account the transition explicitly in this 
comparison with models. Nevertheless, we will explore the effects of a gravity-dependent 
transition between cloudy and cloudless atmospheres in Section 5.4 as this is a topic of 
growing interest. 



The Saumon & Marley (2008) models were computed with what has come to be known 



as a traditional or hot-start initial condition. As discussed in Baraffe et al. (2002), Marley 



et al. (2007a) and Spiegel & Burrows (2012) however, the computed radii of young giant 



planets at ages of 100 Myr and less is highly dependent on the details of the assumed initial 
condition. Even assuming very cold initial conditions, however, computed planetary radii 
never fall below 1 Rj at ages of less than 1 Gyr. Rather than carrying out the model fitting for 
an uncertain set of assumed cold initial conditions, we choose here to employ the traditional 



-14- 



hot-start boundary conditions for the evolution modeling. In this way we avoid unphysical 
very small radii (R < 1 Rj) while adding an additional constraint to the modeling. 



4. Application to HR 8799 Planets 
4.1. Constraints on the HR 8799 System Properties 



A number of the properties of the HR 8799 system as a whole help to constrain the 
properties of the individual planets. Of foremost importance of course is the age of the 
primary star since older ages require greater planetary masses to provide a fixed observed 
luminosity. The massive dust disk found outside of the orbit of the most distant planet, HR 
8799 b, constrains the mass of that planet since a very massive planet would disrupt the 
disk. Finally dynamical models of the planetary orbits circumscribe the parameter space 
of orbits and masses that are stable over the age of the system. All of these topics have 
been discussed extensively in the literature so here we briefly summarize the current state 



of affairs. A more thorough review can be found in Sudol & Haghighipour (2012). 



Since the discovery of the first three planets, the age of HR 8799 has been debated. As 
summarized initially by the discoverers, most indicators suggest a young age of 30 to 60 Myr 



(Marois et al. 2008). However the typical age metrics are somewhat more in doubt than 



usual because HR 8799 is a A Boo-type star with an unusual atmospheric and uncertain 



internal composition. Moya et al. (2010) review the various estimates of the age of the star 



prior to 2010 and argue that most of the applied metrics, including color and position on 



the HR diagram, are not definitive. Most recently Zuckerman et al. (2011) conclude that 



the Galactic space motion of HR 8799 is very similar to that of the 30 Myr old Columba 
association and suggest that it is a member of that group. They also argue that the B — V 
color of HR 8799 in comparison to Pleiades A stars also supports a young age, although 
the unusual composition hampers such an argument. Perhaps the fairest summary of the 
situation to date would be that most traditional indicators support a young age for the 
primary, but that no single indicator is entirely definitive on its own. 

One indicator that the age could be much greater than usually assumed is discussed by 



Moya et al. (2010). Those authors use the 7 Doradus g-mode pulsations of the star to place 



an independent constraint on the stellar age. Their analysis is dependent upon the rotation 
rate of the star and consequently the unknown inclination angle and thus is also uncertain. 
Nevertheless they find model solutions that match the observed properties of the star in 
which the stellar age can plausibly be in excess of 100 Myr and in some cases as large as 1 
Gyr or more. They state that their analysis is most uncertain for inclination angles in the 
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range of 18 to 36°, which corresponds to the likely inclination supported by observations of 
the surrounding dust belt (see below). Thus stellar seismology provides an intriguing, but 
likewise still uncertain constraint. 

The dust disk encircling the orbits of the HR 8799 planets can in principle provide 
several useful constraints on the planetary masses and orbits. First the inclination of the 



disk affects the computed orbital stability of the companions (Fabrycky & Murray-Clay 2010) 



if we assume the disk is coplanar with the planetary orbits. If the rotation axis of the star 
is perpendicular to the disk, the inclination also has a bearing on the stellar age since the 



seismological analysis in turn depends upon its inclination to our line of sight (Moya et al. 



2010). Hughes et al. (2011 ) discuss a variety of lines of evidence that bear on the inclination, 



i, of the HR 8799 dust disk. While they conclude that inclinations near 20° are most likely, 
the available data cannot rule out a face-on (i = 0°) configuration. Finally an additional 
important constraint on the mass of HR 8799 b could be obtained if it is responsible for 
truncating the inner edge of the dust disk. An inner edge at 150 AU is consistent with 
available data (Su et al. 2009) and this permits HR 8799 b to have a mass as large as 20 Mj 



(Fabrycky & Murray-Clay 2010). It is worth noting, however, that this limit depends upon 



the model-dependent inner edge of the disk and the dynamical simulations. 

Finally dynamical simulations of the planetary orbits constrained by the available astro- 



metric data can provide planetary mass limits. In the most thorough study to date Fabrycky 



& Murray-Clay (2010) found that if planets c and d were in a 2:1 mean- motion resonance 



their masses could be no larger than about 10 Mj. However if there were a double resonance 
in which c, d, and b participated in a "double 2:1" or 1:2:4 resonance (originally identified by 



Gozdziewski & Migaszewski (2009)) then masses as large as 20 Mj are permitted and such 
systems are stable for 160 Myr ( |Fabrycky fc Murray-ClaypHO] ). Such a resonance was found 
to be consistent with the limited baseline of astrometric data. HR 8799 b,c, and d have also 



been identified in an archived HST image taken in 1998 (Lafreniere et al. 2009 Soummer et 



al. 2011). These data continue to allow the possibility of the 1:2:4 mean motion resonance, 
a solution which implies a moderate inclination (i = 28°) for the system. New dynamical 
models that include both this new astrometric data and the innermost e planet are now 



required to fully evaluate the system's stability. Sudol & Haghighipour (2012) studied such 



a system with masses of 7, 10, 10, and 10 Mj. They generally found system lifetimes shorter 
than 50 Myr for such large masses but at least one system was found to be stable for almost 
160 Myr. 

Taken as a whole the age of the system and the available astrometric data and dynamical 
models are consistent with a relatively young age (30 to 60 Myr) and low masses for the 
planets (below 10 Mj). However the possibility of an older system age, as allowed by the 
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asteroseismology, and higher planet masses, as permitted if the planets are in resonance and 
by the dust disk dynamics, cannot be fully ruled out. Given this background we now consider 
the planetary atmosphere models. 



4.2. Data Sources 



The available photometric data for each planet is summarized in Table 2 and shown 
on Figures [SHS} In addition for planet b we employ H and K band spectra as tabulated in 



Barman et al. (2011a). We do not include the narrow band photometry of Barman et al. 



(2011a) since this dataset has been superseded by the spectroscopy. We also do not include 



very recent 3.3-/im photometry from Skemer et al. (2012) which became available after the 



submission of this manuscript although we do plot the point in Figures |4]-[6j Below we sum- 
marize the sources of the photometry used in the fitting. With the exception of the Subaru 
z-band which sits in an atmospheric window, we included an atmospheric transmission curve 
when computing the synthetic magnitudes of the model spectra. The transmission curve 
was generated with ATRAN (Lord 1992) at an airmass of 1 with a precipitable water vapor 
content of 2 mm. 



4-2.1. Subaru-z band 



The Subaru- z-band photometry is from Currie et al. (2011) and was obtained with the 



Infrared Camera and Spectrograph (IRCS; Tokunaga et al. (1998)) on the Subaru Telescope 



The filter profile was kindly provided by Tae-Soo Pyo. No atmospheric absorption was 
included because the filter sits in a window that is nearly perfectly transparent. 



4.2.2. J band 



The J band data were taken from Marois et al. (2008) and Currie et al. (2011). The 



former observations were done with the Near-Infrared Camera (NIRC2) on Keck II which 
uses a Mauna Kea Observatories Near- Infrared (MKO-NIR) J band filter. We used the filter 
transmission profile from Tokunaga et al. (2002). The latter observations were obtained 



with the Infrared Camera and Spectrograph (IRCS; Tokunaga et al. (1998)) on the Subaru 
Telescope which also uses a MKO-NIR J band filter. 
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4-2.3. H and Ks bands 



The if -band and f^-band data were taken fromjMarois et al. (2008). The observations 



were done with the Near-Infrared Camera (NIRC2) on Keck II which uses MKO-NIR filters. 



We used the filter transmission profile from Tokunaga et al. (2002) 



4.2.4. [3.3] band 



The [3.3]-band data was taken from (Currie et al. 2011). The observations were done 
with the Clio camera at the MMT Telescope (Freed et al. 2004 Sivanandam et al. 2006). 
The filter is non standard and has a central wavelength of 3.3 /mi, and half-power points of 
3.10 and 3.5 /im. The filter transmission profile was provided by Phil Hinz. 



4.2.5. V band 



The //-band data was taken from Currie et al. (2011). The filter is the V filter in the 
MKO-NIR system so we used the filter transmission profile from Tokunaga et al. (2002). 



4.2.6. M'-band 



The M-band photometry of Galicher et al. (2011) was obtained using the Near- Infrared 



Camera (NIRC2) on Keck II. This filter profile is the same as the M' band of the MKO-NIR 



system. We therefore used the filter transmission profile from Tokunaga et al. (2002) 



4.3. Fitting Method 

In order to determine the atmospheric properties of the HR 8799 planets, we compared 
the observed photometry to synthetic spectra generated from our model atmospheres. We 
used a grid of solar metallicity models with the following parameters: T e g = 600-1300 K in 
steps of 50 K, log g = 3.5-5.5 in steps of 0.25 dex, / so d = 1,2, and eddy mixing coefficient 
K zz = 0, 10 4 cm 2 s _1 . We identify the best fitting model and estimate the atmospheric 
parameters of the planets following the technique described in Cushing et al. (2012, in prep). 
In brief, we use Bayes' theorem to derive the joint posterior probability distribution of the 
atmospheric parameters given the data P(T e g, logg, / se d, K zz \f), where f represents a vector 
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of the flux density values (or upper limits) in each of the bandpasses. Since the posterior 
distribution is only known to within a multiplicative constant, the practical outcome is a list 
of models ranked by their relative probabilities. 

Estimates and uncertainties for each of the atmospheric parameters can also be derived 
by first marginalizing over the other parameters and then computing the mean and standard 
deviation of the resulting distribution. For example, the posterior distribution of T e g is given 
by, 

P(T eS \f) = J P(T eff , log g, f S ed,K zz \f) d logg df sed dK zz 

Since (T eff ,logg) values can be mapped directly to (M, R, L ho {) values using evolutionary 
models, we can also construct marginalized distribution for M, R, and L^ \. Figure [7] shows 
the resulting distribution of T c g, logg, M, and Lboi for each planet and indicates the formal 
solution for these parameters and and their associated uncertainties. 

Finally note that we chose to use a Bayesian formalism rather than the more common 
approach of minimizing \ 2 because 1) we can marginalize over model parameters such as the 
distance and radii of the brown dwarfs, and 2) we can incorporate upper limits using the 



formalism described in Isobe et al. (1986). 



4.4. Results of Model Fitting 

In this section we discuss the individual best fits to each planet. Figures [4] [6] display the 
model fits to the observed spectra and photometry. Each panel of Figure [8] shows contours, 
denoting integrated probabilities of 68, 95, and 99%, in the logg — T c s plane. In these 
figures evolution tracks for planets and brown dwarfs of various masses are shown. The 
objects evolve from right to left across the figures as they cool over time. Isochrones for a 
few ages are shown; the kinks arise from deuterium burning. In some cases at a fixed age 
a given T e s can correspond to three different possible masses (e.g., a 1150 K object at 160 
Myr). Also shown are contours of constant Lboi- Note that the isochrones are derived from 
the conventional hot-start giant planet evolution calculation. A different choice of initial 
conditions would result in different isochrones. 

The best fitting parameters are also shown in Figure [7] as histograms of probability 
distribution for T e g, logg, M and L. For logg and T e g the histograms are projections of the 
contours shown in Figure [8] onto these two orthogonal axes. The mean of the fit and the 
size of the standard deviation is indicated in each panel and also illustrated by the solid and 
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dashed vertical lines. The third and fourth columns of Figure [7] depict the same information 
but for the mass and luminosity corresponding to each (T e ff, logg) pair, as computed by the 
evolution model. 

We discuss each set of fits for each planet in turn below. 



4.4.I. HR 8799b 



HR 8799b is the only one of the three planets considered here for which there is spec- 
troscopic data and our results are sensitive to whether or not this data is included in our fit. 
Contours which show the locus of the best fitting models for the photometry are shown in 
the left-hand panel of Figure [8j When only the photometric data is fit high masses around 
~ 26 Mj are favored. The photometry-only fit finds T e g = 1000 K and / se d = 2 while a fit 
to both the spectroscopy and the photometry results in T eS = 750 K and / se d = 2 with a 
mass of ~ 3Mj. We reject the low temperature fit for several reasons: the solution lies at 
the edge of our model grid, such a planet would be very young, and such a cold effective 
temperature is not consistent with the bolometric luminosity of planet b (see §5.2). These 
models are illustrated in the top two panels of Figure [21 



To isolate the effect of the spectroscopy of Barman et al. (2011a) on the preferred fit 



we relaxed the radius and distance constraint on the fitting and found the model that best 
reproduces the shape of the spectra. Somewhat surprisingly this is a cold, very low gravity 
and very cloudy model (T cS = 600 K, log g = 3.5 and / se d = 1). With a standard radius such 
a model is again too young and faint and also lies at the edge of the model grid. 

The reason the derived gravity depends so strongly on the H and K spectra is that the 
shape of the emergent flux-and not just the total flux in a given band-contains information 
about the gravity. In particular a "triangular" H band shape serves as an indicator of low 
gravity (see Rice et al. (2011) and Barman et al. (2011a)). This shape results from the 



interplay of a continuum opacity source-either cloud opacity (in a cloudy atmosphere) or 
the collision-induced opacity of molecular hydrogen (when cloud opacity is unimportant)- 
and a sawtooth-shaped water opacity (discussions in the literature generally only highlight 
the latter). At high pressures the continuum hydrogen opacity and/or the cloud opacity 
tends to fill in the opacity trough at the minimum of the water opacity in H band. Since the 
photosphere of lower gravity objects at fixed effective temperature is at lower pressures, the 
H2 and cloud opacity is somewhat less important allowing the angular shape of the water 



opacity to more strongly control the emergent flux (see Figure [9] and Figure 6 of Rice et al. 
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Thus we find that the shape of the H band spectrum is responsible for pulling the 
preferred model fits to low gravity and low effective temperature. Weaker methane bands 
at lower logg in this T e g range also push the fit to lower gravity. The greater number of 
datapoints in the spectra overwhelms the photometric data which is why the contours for 
the best overall fit lie outside of the accepted luminosity range. As we discuss in Section 5.1 
our preferred interpretation is that none of our current models match the true composition, 
mass, and age of this planet. 

The model which best fits the photometry alone in the top panel of Figure [4] fits the 



YJHK and [3.3]-/xm (but not the revised Skemer et al. (2012) [3.3]) photometry to within 



la. The model is too bright at L' and M' . The photometry plus spectrum fit features a 
methane band head at 2.2 /zm that is too prominent, even with log K zz = 4. Both sets of 
solutions, are inconsistent with the accepted age of the the star. The lower mass solution 
would imply very young ages for the planet, well below 30 Myr. Conversely the higher mass 
range implies ages in excess of about 300 Myr. Thus along with the discarded low mass 
fit the photometry-only, higher mass fit is problematical since the mass conflicts with the 
constraints discussed in Section 4.1 



44.2. HR 8799c 

For planet c there is no available spectroscopy and we fit only to the photometry. The 
formal best fitting solution yields T e g = 980 ± 70 K and logg = 4.33 ± 0.28 for a mass of 
15±8Mj. However in both the contour diagram (Figure [8]) and the histogram (Figure [7]) we 
find two islands or clusters of acceptable fits, one at higher gravity and effective temperature, 
and one with lower values for both. The high mass solution lies at masses greater than 20 Mj 
and T e g ~ HOOK. Such models are consistent only with ages around 300 My, well in excess 
of the preferred age range for the primary and the dynamical constraints on the mass. The 
second island of acceptable fits lies at logg ~ 4.25 and T e g ~ 950 K. Figure [5] illustrates the 
spectra for the best fitting model from each case. The lower mass model has logg = 4.25, 
f se d = 2, and logK zz = 4, implying M ~ 10 Mj which is consistent with the dynamical mass 
constraint and represents our preferred solution and is listed in Table 1. The age predicted 
by the evolution of these models is about 160 Myr, consistent with the asteroseismological 
age constraint but not the generally favored range of 30 to 60 Myr. However models with 
modestly lower gravity and slightly smaller masses also fall within the la contours seen in 
Figure [8] do lie within this age range. 

The cooler model fits most of the photometric points to within 2a or better, but varies 
most significantly from the data at [3.3] /im and L', which perhaps imply that despite the 
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disequilibrium chemistry the models have too much methane. The lower gravity solutions 
differ from the high gravity ones most prominently in the red side of K band (where the cooler 
model has a much more prominent methane band head) and at 3 to 4 /im. By constraining 
the methane band depth in the K band and, to a lesser extent, in the H band, spectroscopy 
has the potential to distinguish between these two cases. The shape of H band (Figure 
[9j) can also serve as a gravity discriminator with a more triangular shape indicating lower 
gravity. 



44.3. HR 8799d 

Because of larger observational error bars, the model fits for the innermost of the three 
planets considered here are the most uncertain. As seen in Figure [8] the best fitting models 
allow masses ranging from 5 to 60 Mj and T e g between 900 and 1200 K. However the very 
best fitting models favor solutions with logg around 4.25 to 4.50 and T eff = 1000 K yielding 
a mass of 10 to 20 Mj. As with planet c such a solution is consistent with the dynamical 
constraint but not the age constraint. Also as with planet c the lower end of this mass range 
offers marginally poorer fits that nevertheless still lie within the la contour and that do 
satisfy the age constraint. The best fitting spectrum is shown in Figure [6j 



5. Discussion 

5.1. Implied Masses and Ages 

To summarize our findings from the previous section, each of the three planets considered 
presents a different challenge to characterize. Some model fits to planets c and d imply im- 
plausibly large masses or ages but other acceptable fits satisfy all of the available constraints. 
Both c and d can be characterized as having masses as low as 7 to 8Mj, T eff = 1000 K, and 
fsed — 2 which implies ages of around 60 Myr, within the most commonly cited age range of 
the primary. Some better fitting models have slightly larger masses (10 Mj) and ages (160 
Myr). This age is greater than the range of ages typically quoted for the primary star of 
30 to 60 Myr, although it is within the range permitted by the asteroseismology. Evolution 
models starting from a cooler initial state than the hot-start models would reach these effec- 
tive temperatures and gravity at a younger age than 160 Myr and be more in accord with 
the usual age range. 

For planet b none of the models are satisfactory. Since we do not allow arbitrary radii 
models to fit the data (with the exception of the lowermost panel in Figure pj, we cannot 
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invoke what we judge to be unphysical radii to produce acceptable fits. The solution which 
best fits the photometry alone has M = 26 Mj, T cff = 1000 K, and / se d = 2, but this 
mass clearly violates the constraints discussed in Section 4.1. A fit to the entire spectral 
and photometric dataset results in M ~ 3Mj, T e g = 750 K, and / se d = 1- However we 
discard this model as discussed in Section 4.4.1. This effective temperature is cooler than 



favored by Barman et al. (2011a) and Currie et al. (2011) but is comparable to that found 



by Madhusudhan et al. (2011). 



The most likely explanation for the difficulty in fitting this object is that one of the 
assumptions of the modeling is incorrect. Barman et al. (2011a) speculate that a super- 



solar atmospheric abundance of heavy elements might explain the departures of the data 
from the models. Indeed all of the atmospheres of solar system giant planets are enhanced 
over solar abundance with a trend that the enhancement is greater at lower masses. For 
example Saturn's atmosphere is enhanced in methane by about a factor of ten while Jupiter 



is only a factor of about three (see Marley et al. (2007b) for a review). The available data 
on exoplanet masses and radii suggest that lower mass planets are more heavily enriched in 



heavy elements than higher mass planets (Miller & Fortney 2011). If the mass of HR 8799b 



is intermediate between our two sets of best fits, for example with a mass near 6 or 7Mj, as 
favored by the discovery paper, and if atmospheric abundance trends are similar in the HR 
8799 system to our own, then it may not be surprising if planet b has different atmospheric 
heavy element abundances than c and d. We will consider non-solar abundance atmosphere 
models in a future paper. The full range of model phase space has certainly not yet been 
explored. 

Overall we find that a consistent solution can be found for planets c and d in which 
both have similar masses and ages. This is essentially the solution favored by the discovery 



paper 


Marois et al. 


2008 


) and is within the 


et al. 


( 


2011 


) and 


Madhusudhan et al. 


(2011) 



studies in our finding that the radii for planets b and c that are fully consistent with that 
expected for their individual masses. Unusual radii are not required. 



5.2. Bolometric Luminosities 



The distance to HR 8799 has been measured as d = 39 ± 1.0 pc (van Leeuwen 2007) 



and thus the bolometric luminosity of each planet can be computed from the observed 



photometry. In the discovery paper, Marois et al. (2008) compare the photometry available 



at that time to models and brown dwarf spectra and report the now commonly cited results 
logL bol /L = —5.1 ± 0.1 for planet b and — 4.7 ± 0.1 for c and d. 
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Since the work of Marois et al. (2008), the photometry of the three planets has been 
expanded to cover the SED from ~ 1-4.8 /xm. This better constrains Lboi as ~80% of the 
flux is emitted at these wavelengths. In principle, the bolometric luminosity can be obtained 
by fitting synthetic photometry to the observations, with a scaling factor chosen to minimize 
the residuals. The integrated scaled flux of the model and the known distance gives Lboi 



(Marois et al. 2008). The fitted model thus provides an effective bolometric correction to 



the photometry by approximating the flux between the photometric bands. The scaling 
factor corresponds to (R/d) 2 , where R is the radius of the planet. The optimized scaling 
thus corresponds to an optimization of the radius independent of the physical radius of the 
planet. As is well known, this results in radii for the HR 8799 planets that are considerably 
smaller than be accounted for with the evolution models (Section 1.1). The approach can 
also lead to unrealistic bolometric corrections if the fitted T e g deviates too far from the actual 
value. 

To circumvent this difficulty, here we determine Lboi by using the radius obtained from 
our evolution sequences, which is consistent with our approach to fit the photometry. Of 
course such theoretical radii have their own uncertainty, including a dependence at young 



ages - particularly below 100 Myr - on the initial conditions (Baraffe et al. 2002; Marley et 



al. 2007a Spiegel & Burrows 2012). We neglect the dependence on initial conditions since 



planets forming in the 'cold-start' calculation of Marley et al. (2007a) never get as warm 



or as bright as the HR 8799 planets. Intermediate cases, such as explored by |Spiegel EE 



Burrows (2012) are possible, but we set those aside for now. Our approach, however, does 



eliminate unphysical solutions by constraining the radius to reasonable values (in excess of 
lRj). Thus, for each fitted model (T e fj , log g) we obtain a Lboi from the radius R(T e ff, log g) 
obtained with the evolution^] The resulting probability distributions of L bo i for each planet 
are shown in Fig. [7| along with the mean value and dispersion of each distribution. 

Our fits are based on a model grid with spacing of 50 K and 0.25 dex in T e g and logg, 
respectively, which introduces an additional uncertainty intrinsic to the fitting procedure of 
about half a grid spacing, or ±25 K and ±0.13 dex. We derive the corresponding uncertainty 
in Lboi as follows. The bolometric luminosity is given by 



^bol 



AnGM Q aTfff / M 



off 



g vM 

where the symbols have their usual meaning. From the cloudy evolution of Saumon fc Marley) 



(2008), we find an approximate relation for M(T e g, log g) in the range of T e g and mass of 



2 With "perfect" atmosphere and evolution models the two methods would give identical results. 
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interest: 



log 



M 

Mo 



0.746 log# + 



T, 



<-ir 



5090 



5.35, 



where T e g is in K and g in cms 2 . Thus, 



log L bol = 4 log T efr + 



cff 



5090 



0.254 log g + A, 



where A is a constant. With the grid spacing uncertainties given above, we find A logL bol = 
±0.054, which we round up to 0.06. Combining quadratically this uncertainty with the 
dispersion in Lboi found in our fits (Fig. [7]), we find the luminosity for planet b to be 
log L hol /L e = -4.95 ± 0.06, -4.90 ± 0.10 for planet cR and -4.80 ± 0.09 for planet d. 



These values are consistent with those reported by Marois et al. (2008) although they are 
0.1 dex brighter for planet b, 0.2 dex fainter for planet c, and 0.1 dex fainter for planet d. 
The quoted uncertainties are lower limits of course, since they do not account for obvious 
systematic errors in the models (Figures El-p]). 



5.3. Cloud Properties 



Although there is a dispersion in the best fitting logg and T e ff, essentially all of the 
acceptable fits require a cloud sedimentation efficiency of / se d = 2. As shown in Figure 1 
this value is typical of the best fitting parameters for most field L dwarfs we have previously 



studied (Cushing et al. 2008 Stephens et al. 2009). The persistence of clouds to lower 



effective temperatures at low gravity is also apparent from this figure. By 1000 K most field 
dwarfs with log g > 5 have already progressed to / se( j > 4 whereas clouds persist much more 
commonly among lower gravity objects down to 1000 K. By very cool effective temperatures, 
however, the silicate and iron clouds have certainly departed from view as demonstrated by 
the one logg = 4, T eff ~ 500 K object (ULAS J133553.45+113005.2, ( |Burningham et al. 
2008] |Leggett et al.|[2509| )). 



As Figure 1 attests, the cloud in planets b, c, and d are unusual not so much for their 
global characteristics (the same cloud model that describes L dwarf clouds fits them as well), 
but rather for their persistence. At / se d = 2 there are three field objects with T e fj < 1200 K. 



These objects are 2MASS 0825196+211552 (Kirkpatrick et al. 2000), SDSS 085758+570851 



(Geballe et al. 2002), and SDSS J151643.01+305344.4 (Chiu et al. (2006); hereafter SDSS 



1516+30). Their infrared spectral types are L6, L8, and TO. 5 and the first two are both 
redder in J — K than is typical for those spectral types (Stephens et al. 2009). 



'Note that the dispersion for planet c is non-Gaussian (Fi 
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Figure 3 compares some of the model silicate cloud properties of a low gravity planet 
with models for field L6 and TO. 5 objects. As expected from the discussion in Section 2.2, the 
lower gravity model is marked by a larger particle size than the higher gravity models, and 
the column optical depth of the silicate cloud in all three objects ends up being very similar. 
More importantly the range of cloud optical depths that lie in the near-infrared photosphere 
are similar for all three objects. Thus a low gravity (log g = 3.5) object with T e g- = 1000 K 
ends up with cloud opacity that is very similar to a high gravity (log g = 5.5) object with 



T e s = 1200 K and consequently similar spectra and colors. Indeed Barman et al. (2011a) 
has already noted the similarity of SDSS 1516+30 to HR 8799b. This congruence between 
lower gravity and higher gravity models led to the initial surprise that the apparently cool 
planets seem to have clouds reminiscent of higher gravity-and warmer-L dwarfs. 

The relative contribution of clouds to the opacity in individual photometric bands is 



depicted in Figure [TOj This figure presents contribution functions for the J, H, K, L', 
and M' bands for six different combinations of gravity, effective temperature, and cloud 
treatment. The contribution functions illustrate the fractional contribution to the emergent 
flux as a function of pressure in the atmosphere. In a cloud-free, T eS = 1000 K, logg = 5.0 



atmosphere (left panel, Figure 10 1) the V flux emerges predominantly near P = 0.6 bar 
while the J-band flux emerges from near 8 bar. The contribution functions do not account 
for the effect of cloud opacity, but rather show for each case where the flux would emerge 
from for that particular model if there were no clouds. 



The center two panels of Figure 10 1 and b illustrate the vertical location of the cloud 
layers for both / se d = 1 and 2. The / se d = 2 clouds are thinner and the cloud base is deeper 
since these less cloudy atmospheres are cooler than the f se & = 1 case, as seen in the right 
hand panels. If the cloud deck lies above or overlaps the plotted contribution function of 
a given band then the emergent flux in that band will be strongly affected by the presence 
of the cloud. The figure makes clear that regardless of gravity thicker clouds impact more 
of the emergent spectra than thinner clouds. Clouds described by / se d = 2 strongly impact 
J, H, and K bands, but are less important at L', and M'. We conclude that at least for 
the effective temperature range inhabited by HR 8799 b, c, and d that clouds are most 
strongly impacting the observed spectra at wavelengths shorter than about 2.5 /xm while the 
longer wavelength flux is primarily emerging from above the cloud tops. Figures such as 
this illustrate the value multi-band photometry has in both constraining not only the total 
emergent flux, but also the vertical structure of the clouds. 
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5.4. Evolution with a gravity-dependent L to T transition 



The growing evidence that the cloudy to cloudless transition in field brown dwarfs 
depends on gravity (§2.1) is complemented by the published analyzes of the HR 8799 planets 
(including the present work) which all indicate that their atmospheres are cloudy and that 
they have T e g well below the estimated ~ 1400 K limit of the L dwarf sequence. Thus, it 
appears that the atmospheres of lower gravity dwarfs and of imaged exoplanets retain their 
clouds to lower T e g, which is supported by simple cloud model arguments (§2.2). As we have 
argued, this is the simplest interpretation of the fact that the HR 8799 planets have T e g 
typical of cloudless T dwarfs but have evidently cloudy atmospheres. How is the evolution 
of brown dwarfs across the transition from cloudy to clear atmosphere affected? 

The atmosphere of a brown dwarf largely controls its evolution because it acts as a 
surface boundary condition for the interior. A more opaque atmosphere (more clouds, or 
higher metallicity, for instance) slows the escape of radiation and increases the cooling time 



of the interior. Saumon & Marley (2008) looked at the evolution of brown dwarfs across the 
transition by assuming that the atmosphere was cloudy (/ se d = 2) down to T e g = 1400 K, and 
clear below 1200 K, with an linear interpolation of the atmospheric boundary condition in 
the transition regime. This effectively corresponds to increasing the sedimentation efficiency 
across the transition, one of the proposed explanations for the cloud clearing (§2.1). By 
converting the evolution sequences to magnitudes using synthetic spectra (f se d = 1 for cloudy 
atmospheres, and / se d = 4 for "clear" atmosphere^ a good match to the near-infrared color 
magnitude diagrams of field dwarfs was found from the cloudless late M dwarfs, along the 
cloudy L dwarf sequence, across the L/T transition and down to late T dwarfs. 

We now extend this toy model to include a gravity-dependent range of T e s for the 
transition to explore the consequences, at the semi-quantitative level, on the cooling tracks 
of brown dwarfs and exoplanets. In view of the success obtained for field dwarfs (of relatively 



high gravity) with the Saumon & Marley (2008) toy model, and the requirement that the 
lower gravity HR 8799 planets be cloudy at T e g ~ 1000 K, we define the transition region 
to be T cS = 1400 to 1200 K at \ogg = 5.3 (cgs) and 900 to 800 K at \ogg = A with a linear 



interpolation in between (Fig. 11). The cloudy boundary condition above the transition is 
based on our / sed = 2 atmosphere models, and our cloudless models below the transition, as 
in Saumon & Marley (2008). Synthetic magnitudes are generated from the cooling tracks 



using our new / se d = 1 and cloudless atmosphere models (Saumon et al. 2012). 



4 These are not fully consistent with the values used for the evolution, but the effect on the evolution of 
this small difference in / se d is small. 



-27- 



The resulting cooling tracks of two low-mass objects of 5 and 20 Mj are shown in Fig. 
12 where the same calculation, but based on a fixed T e g transition (Fig. 11) is also displayed 



for comparison. It is immediately apparent that these low-mass objects, which retain their 
clouds to lower T e s (~ 850 K for 5Mj and ~ 1050 K for 20 Mj) with the prescribed gravity- 
dependent transition evolve along the L dwarf sequence longer and reach the region of the 
color-magnitude diagram occupied by the HR 8799 planets before they turn to blue J — K 
colors as the cloud clears. Also remarkable is that in the transition region where the J — K 
color changes from ~ 2 to ~ 0, the low mass object is fainter in K than the higher mass 
object, the reverse of the situation for a transition that is independent of T e g. This effect 
persists up to a cross over mass of ~ 60 Mj above which the trend reverses (Fig. fTTJ. This 
implies that low mass objects that are in the transition region should appear below (i.e. be 
dimmer) the field T0-T4 dwarfs, perhaps by up to 1-2 magnitudes. We note that the pile 
up of objects in the transition region reported in Saumon & Marley (2008) still occurs in 
this new calculation but it is more spread out in T e g, as would be expected from the broader 
span of the transition in T e g (Fig. 11). 

We emphasize that this evolution calculation is a toy model that has been loosely ad- 
justed to account for limited observational constraints. It reveals trends but is not quantita- 
tively reliable. In particular, we have had to use f se d = 1 to match the near infrared colors 
of the HR 8799 planets while our best fits give / se d = 2 for all three planets. This reflects 
the fact that the models give different best-fit parameters when applied to a subset of the 
data, a well-known difficulty with current models (Cushing et al. 2008 Patience et al.|[201~2). 



5.5. Mixing 



Given the discussion in Section 1.3 regarding the prevalence of atmospheric mixing 
resulting in departures from chemical equilibrium in solar system giants and brown dwarfs, 
it is not surprising that mixing is also important in warm exoplanet atmospheres as well. 



Barman et al. (2011a) discuss the influence of non-equilibrium chemistry at low gravity and 



find that the CO/CH4 ratio can become much larger than 1 in the regimes inhabited by the 



HR 8799 planets. Also Barman et al. (2011b) found non-equlibrium chemistry was likely 
important in 2M1207b. 

We find that all of the best fitting models for each planet, b, c, and d, include non- 
equilibrium chemistry. Within our limited grid with K zz = and 10 4 cm 2 s~ 1 , the latter 
choice was strongly preferred in all cases providing yet another indication of the importance 
of chemical mixing in substellar atmospheres. This also suggests that a fuller range of models 
with a greater variety of eddy mixing strengths should be considered in future studies to 
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better constrain this parameter. 



5.6. Mechanism for Gravity Dependent Transition 

In Section 2.2 we demonstrated that the effect of a given cloud layer, all else being equal, 
is greater in a lower mass extrasolar giant planet than in a more massive brown dwarf of 
the same effective temperature. If we add effective temperature as a variable then we find 
that a cooler low mass object can have clouds comparable to a warmer high mass object. 
Such a congruence is empirically demonstrated by the similar spectra of SDSS 1516+30 and 



HR 8799 b (as originally noted by (Barman et al. 2011a)). The former is a ~ 70 Mj, 1200 



K field L dwarf while the latter is plausibly a few Jupiter mass, 1000 K young gas giant 
planet (although the modeling discussed here does not select this solution). Likewise in 
Section 5.4 our simple evolution calculation with a gravity- dependent L to T type transition 
temperature illustrates that the location of young objects on the color magnitude diagram 
can be understood if clouds remain to lower effective temperatures at lower gravity. The fact 
that such behavior is dependent upon gravity is not in itself surprising as a lower gravity 
would be expected to alter its behavior. However the specific question remains, what is 
the specific mechanism that results in lower mass objects making the L to T type spectral 
transition at lower effective temperatures than higher mass objects? In this section we offer 
some speculation while recognizing that a serious analysis is beyond the scope of this paper. 

A possible contributing factor might be found in the relative positions of the convec- 



tion zone and the photosphere as a function of gravity (a point also raised in Barman et 



al. 


(2011a|b 


) and 


Rice et al. 


(2011 



functions for different bandpasses are shown for two different gravities. At T c g = 1000 K 
for moderately cloudy (/ se d = 2) atmospheres the convection zone, regardless of gravity, 
penetrates into the cloud layers that control the J and H band fluxes. For cloudless atmo- 
spheres, however, the convection zone for the high gravity case is quite deep (P > 20 bar), 



well below even the region probed by the J band (Figure 10 1). At lower gravity, however, 
the convection zone penetrates higher into the atmosphere to much lower pressure, over- 
lapping the J band contribution function (Figure [To|d). If we imagine that a given patch 
of atmosphere begins to clear, perhaps because of more efficient local sedimentation, in the 
high gravity case the removal of cloud opacity leads the atmosphere to become radiative and 
more quiescent, favoring particle sedimentation relative to convective mixing and enlarging 
what began as a localized clearing. At low gravity however the removal of cloud opacity does 
not as dramatically push the atmosphere to a quiescent state. Thus convection continues 
to loft cloud particles and the local clearing fills back in. Only when the clear atmosphere 
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convection zone lies very deep do the clouds dissipate. Since low gravity atmospheres are 
more opaque than high gravity ones this process of the growth of clearings begins at lower 
effective temperature at lower gravity. 

Another possibility is that detached convection zones play a role in hastening the L to T 
transition. Within some effective temperature ranges there are two atmospheric convection 
zones, one deeply seated and a detached zone that is separated by a small radiative zone. 



This can be seen in the / se d = 1 temperature profiles in Figure 10 Burrows et al. (2006) and 



Witte et al. (2011) have speculated that the interplay of dynamical and cloud microphysics 
effects that may occur when the intermediate radiative zone forms or departs may play a 
role in the transition. Perhaps at some effective temperature threshold particles forming in 
the upper convective zone grow large enough that they fall all the way through the cloud 
base and the intermediate radiative zone before they completely evaporate. Depending 
on the efficiency of mixing in the radiative zone this could result in a net transport and 



sequestration of condensate away from the near-infrared photosphere. Witte et al. (2011) 



discuss a similar idea of the convection "fanning" the fall of particles away from the upper 



zone. As seen in Figure [10} however, for both the / se( j = 2 and the cloudless case there is 
only one convection zone, so the potential for multilayered convection is less compelling in 
this case. Nevertheless such mechanisms require more sophisticated modeling to ascertain 
how they might be affected by gravity and effective temperature. 

Arguments such as these that are based upon ID radiative convective models only 
scratch the surface of the underlying complex dynamical problem. For example |Freytag"et 



al. (2010) performed two-dimensional radiation hydrodynamic simulations of brown dwarf 
atmospheres to study the effects of clouds on atmospheric convection. They found that 
atmospheric mixing driven by cloud opacity launches gravity waves that in turn play a 
role in maintaining the cloud structure. The Freytag et al. study considered a domain a 
few hundred kilometers wide by about 100 km deep and only investigated a single gravity 
(log g = 5) so how such effects might vary with gravity is not yet known. Furthermore how 
the local clouds might interact with the very large scale planetary circulation has not been 
explored. Perhaps clouds form holes or otherwise dissipate only when most of the cloud 
optical depth lies deeply enough to be strongly influenced by global atmospheric circulation. 
Large scale global dynamical simulations that capture the relevant physics of particle and 
energy vertical and horizontal transport are likely required to fully describe the L to T 
transition mechanism. 
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5.7. Future 



Our experience in fitting the spectra of planet b in particular points to the importance 
of spectra in the analysis. Adding the H and K band spectra to the analysis results in much 
lower preferred masses than fitting photometric data alone. Thus we expect that additional 
spectral data will further inform future model fits. 

As noted in Section 2.1 one hypothesis for the nature of the L to T transition is that 
it involves partial clearing of the assumed global cloud cover. It is possible that models 
which include partial cloudiness may better describe the observed flux and |Currie et al. 



(2011) have explored this possibility. Given the limited data available today we feel the 



addition of another free model parameter is premature and in any event we have found that 
brown dwarfs with partial cloud cover have an overall near-infrared spectrum that resembles 



a homogeneous dwarf with a thinner, homogenous global cloud (Marley et al. 2010). 



Another method for characterizing these planets and probing atmospheric condensate 



opacity in self-luminous planets is by polarization (Marley & Sengupta 2011; de Kok et al. 



2011). Marley & Sengupta (2011) found that rapidly rotating, homogenously cloud-covered 



planets may be sufficiently distorted to show polarization fractions of a few percent if they are 



relatively low mass, de Kok et al. (2011 ) found that even when partial cloudiness is considered 



much larger polarization fractions are unlikely. However if this level of polarization could 
be measured in one of the HR 8799 planets this would confirm the presence of clouds and 
also place an upper limit on the planetary mass. Objects in this effective temperature range 
(near 1000 K) and with logg > 4 are predicted to exhibit polarization well below 0.2%. Both 
SPHERE and GPI have polarization imaging modes, but it is not clear if they would have 
sufficient sensitivity to place useful upper limits on the HR 8799 system. 



6. Conclusions 

We have explored the physical properties of three of the planets orbiting HR 8799 by 
fitting our standard model spectra to the available photometry and spectroscopy. Unlike 
some previous studies we have required that models with a given \ogg and T e g have a 
corresponding radius that is calculated from a consistent set of evolution models. While 
the radii of the planets are not variables, we do include two other free parameters: the 
cloud sedimentation efficiency / se d and the minimum value of the atmospheric eddy mixing 
coefficient K zz . 

In agreement with all previous studies we find that the atmospheres of all three planets 
are cloudy, which runs counter to the expectation of conventional wisdom given their relative 
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low effective temperature. However as we argue in Sections 2.1 and 2.2, finding clouds to 
be present at lower effective temperatures in lower gravity objects is fully consistent with 
trends already recognized among field L and T dwarfs and from basic atmospheric theory. We 
uniformly find that the best fitting value of the sedimentation efficiency / se d is, in essentially 
all cases, 2, which is typical of the value seen in pre-L/T transition field L dwarfs (Fig. 1) 



(Cushing et al. 2008; Stephens et al. 2009). In agreement with Barman et al. (2011a) we 



thus find that the clouds in these objects are neither "radically enhanced" (Bowler et al. 



2010) nor representative of a "new class" (Madhusudhan et al. 2011) of atmospheres. 



As have some previous authors (Barman et al. 2011a[b ) we find that eddy mixing in 
nominally stable atmospheric layers is an important process for altering the chemical com- 
position of all three planets. While we have not carried out a comprehensive survey of 
non-equilibrium models, we find that values of the eddy mixing coefficient near log K zz ~ 4 
generally fit the available data better than models that neglect mixing. Such values are typ- 
ical of those found for field L and T dwarfs (e.g., Stephens et al.|[2009 ) and the stratospheres 



of solar system giant planets (e.g., see the detailed discussion for Neptune in Bishop et al. 



(1995)). 



The best fitting values for the primary model parameters log g and T e fj are less secure. 



For HR 8799 b the inclusion of the H and K band spectra of Barman et al. (2011a) drive 
our fits to low masses of ~ 3Mj and effective temperatures, a solution which we discard 
as discussed in Section 4.1.1. The photometry alone favors much higher masses, ~ 25 Mj 
that are apparently ruled out by dynamical considerations. Thus we find no plausible model 
that fits all of the accepted constraints. Fits for the planets c and d likewise generally favor 
higher masses, although there are some solutions that are consistent with masses near or 
below ~ 10 Mj with ages consistent with the available constraints. For all three planets the 
photometry predicted by the best fitting model is generally consistent with the observed 
data within 1 to 2 standard deviations. We stress that all of these fits have radii that are 
appropriate for the stated effective temperature and gravity. 

In conclusion the modeling approach that has successfully reproduced the spectra of field 
L and T dwarfs seems to also be fully applicable to the directly imaged planets. Nevertheless 
a larger range of model parameters, including non-solar metallicity, must be explored in order 
to fully characterize these objects as well as the planets yet to be discovered by the upcoming 
GPI, SPHERE, and other coronagraphs. 
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Table 1. Summary of Derived Planet Properties 



Planet 


Ref. 1 


M (M Jup ) 


log g 


T eff (K) 


R(R 3 ) 


age 


(Myr) 


log L bo i/ L & 


b 2 


Blla 


0.1-3.3 


3.5 ±0.5 


1100 ± 100 


0.63 - 0.92 


30 


-300 


-5.1 ±0.1 3 




Cll 


5-15 


4-4.5 


800 - 1000 




30 


-300 






Gil 


1.8 


4 


1100 


0.69 










Mil 


2-12 


3.5-4.3 


750 - 850 




10 


- 150 






M12 4 


26 


4.75 


1000 


1.11 


360 


-4.95 ±0.06 


c 


Cll 


7-17.5 


4-4.5 


1000 - 1200 




30 


-300 


-4.7±0.1 3 




Gil 


1.1 


3.5 


1200 


0.97 










Mil 


7-13 


4-4.3 


950 - 1025 




30 


- 100 






M12 


8-11 


4.1 ±0.1 


950 ± 60 


1.32 - 1.39 


40 


- 100 


-4.90 ±0.10 


d 


Cll 


5 - 17.5 


3.75-4.5 


1000 - 1200 




30 


-300 


-4.7±0.1 3 




Gil 


6 


4.0 


1100 


1.25 










Mil 


3-11 


3.5-4.2 


850 - 1000 




10 


- 70 






M12 


8-11 


4.1 ±0.1 


1000 ± 75 


1.33 - 1.41 


30 


- 100 


-4.80 ±0.09 



Currie et al. 


( 


2011 


); Gil = 



Galicher et al. (2011); 



Mll=Madhusudhan et al. (2011); M12=this work. 



2 Parameters derived by Bowler et al. (2010) are not listed because of very large scatter de- 
pending upon various assumptions. 



Luminosity from Marois et al. (2008). 

4 For b this is the formal best fit single model to the photometry alone, for c and d these are the 
preferred solution ranges as discussed in the text. The b fit is incompatible with the generally 
accepted constraints as discussed in the text. Formal solutions are shown in Figure [7| 
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Table 2. Photometric Data for the HR 8799 Planets 



Planet 


TTi 1 

Band 


A 1 TV T 

Abs. Mag. 


Ret. 1 


b 


Subaru-z 


18.24 ± 0.29 


Cll 






1fi ^2 + 14 


CI 1 

V X X 




1 1 


ItiO I _C U.l 1 


M08 




K 


Iti.Uu _l_ W.W(J 


M08 




[3 31 

[O.OJ 


1 3 96 + 98 


V. X X 




TJ 


1 2 68 ± 1 2 


Cll 

V... X X 




1V1 


~\ o D7 _i_ n on 
lo.U f ± U.oU 


pi 1 


c 


Subaru- z 


> 16.48 


Cll 




J 


14.65 ±0.17 


M08 




H 


13.93 ±0.17 


M08 




K s 


13.13 ±0.08 


M08 




[3.3] 


12.64 ±0.20 


Cll 




V 


11.83 ±0.07 


Cll 




M' 


12.05 ±0.14 


Gil 


d 


Subaru-z 


> 15.03 


Cll 




J 


15.26 ±0.43 


M08 




H 


13.86 ±0.22 


M08 




K s 


13.11 ±0.12 


M08 




[3.3] 


> 11.63 


Cll 




M' 


11.67 ±0.35 


Gil 



^ll HCurrie et al.| ( |2011[ ) 
M08 HMaro"is et al.| ( J2008| ) 
Gll=|Currie eFaL] pHl) 
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Fig. 1. — Model parameters / se d and T eff as derived by various applications of Marley & 
Saumon atmosphere and evolution models. Size of dot reflects derived logg(cms~ 2 ) and 
'nc' denotes cloudless models (note that 'nc', which corresponds to / se( j — > oo, is arbitrarily 
plotted at /sod = 5). Points which would otherwise overlap are slightly offset vertically and 
the T e ff values decrease to the right to suggest evolution in time. The points for HR 8799 
c and d from the analysis here are labeled with planet designator. Remaining points are 



from 


Geballe et al. 


(2001 


); 


Mainzer et al. 


(2007) 


; Leggett et al. 


(2007a 


2008 


); 


Geballe et 


al. 





>009 


); 


Leggett et al. 


(2009 


); 


Stephens et al. i 


2009 


); 


Mainzer et al. 





2011) although fits 



to unresolved binaries and objects with very poorly constrained properties (e.g., Gl 229 B 
with logg uncertain by a full dex) are excluded. SDSS 1516+30 is denoted by '1516'. The 
cross denotes size of the typical uncertainties in the model fits which are usually ±100 K in 
effective temperature, ±0.25 dex in log g, and ±0.5 in / sed , although the uncertainty analysis 
is not uniform across the various sources. 
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Fig. 2. — Model atmosphere temperature-pressure profiles for cloudy brown dwarfs and 
planets assuming / sec j = 2 (Ackerman & Marley 2001). Each profile is labeled with log g 
and T e ff of the model. The condensation curve for forsterite is shown with a dotted line. 
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Fig. 3. — Silicate cloud properties as computed by the Ackerman & Marley (2001) cloud 
model for three models. From left to right the the best-fitting models Stephens et al. (2009) 
for 2MASS 0825+21 and SDSS 1516+30 are shown along with a profile for a young, cloudy, 
three Jupiter mass planet. Labels underneath each object name denote model T e fj(K) / 
logg(cgs) / f se d- Dashed curves show the effective radius, r e g of the particles on the top 
axis. The column optical depth as measured from the top of the atmosphere is shown by 
the solid lines and the scale on the bottom axis. Thicker lines denote the region of the cloud 
which lies within the A = 1 to 6 /zm photosphere. Other modeled clouds are not shown for 
clarity. 
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Fig. 4. — Observed (black) and model (red, green, purple) photometry and spectra (see 
Table 1 and Barman et al. (2011a)) for HR 8799b. Models are identified in the upper left 
hand corner of each panel by T cS / \ogg (cgs)/ fsed/K zz . The top panel shows the model that 
best fits the photometry alone while the middle panel shows the solution that best fits both 
the photometry (excluding H and K bands) and spectroscopy simultaneously. Model fluxes 
and photometry have been computed for radii specific to the T e g and log g of the atmosphere 



model at a distance of 39.4 pc as observed from Earth. The [3.3] /im photometry of Skemer 



et al. (2012) is shown as a blue star and is not included in the fits but rather is shown 



for comparison purposes only. The lower panel shows the model that best fits the H and 
.fT-band spectrum alone. However in contrast to the top two panels where the absolute flux 
level of the models are set by the model radii and known distance to HR 8799, the absolute 
flux level of the model in the lower panel is determined by minimizing \ 2 between the models 
and data. 
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Fig. 5. — The two best fitting model spectra for HR 8799 c. Observed photometry (see 
Table 2) is shown in black, high and low gravity solutions in green and red, respectively. 
The two solutions correspond to the centers of the two best fitting islands in the contour plot 
shown in the middle panel of Figure |8j Models are identified in the upper left hand corner 
by T c ff/ log g (cgs)// se( j /K zz . The [3.3] fim photometry of Skemer et al. (2012) is shown as a 
blue star and is not included in the fits but rather is shown for comparison purposes only. 
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Fig. 6. — The best fitting model for HR 8799 d. Observed photometry (see Table 1) is shown 
in black; model photometry is indicated by the red dots. Model is identified in the upper left 
hand corner by T e g/\ogg (cgs)/ f se d/K zz . The 3.3-/im photometry of Skemer et al. (2012) 
is shown as a blue star and is not included in the fits but rather is shown for comparison 
purposes only. 
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Fig. 7. — Histograms depicting the probability density distributions of the various model 
parameters to planets HR 8799 b, c, and d. For planet b only the results for the fitting of 
the photometry are shown. The T e g and logg histograms can be thought of as the projection 
of the contours shown in Figure [8] onto these two orthogonal axes. In each case the mean of 
the fit and the standard deviation are indicated by /i and a, respectively. These quantities 
are in turn illustrated by the solid and dashed vertical lines. For the parameters for planet 
b, only a single model is identified so no standard deviation is given. The third and fourth 
columns of histograms depict the same information as the first two, but for the mass and 
luminosity corresponding to each (T e g, log g) pair, as computed by the evolution model. 
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Fig. 8. — Contours illustrate domain of best-fitting models on the log g — T e g plane. For 
each planet three contours are shown which correspond to integrated probabilities of 68, 



95, and 99% (red, thick to thin contours). Evolution tracks from Saumon et al. (2007) are 
shown as labeled black curves; planets evolve from right to left with time across the diagram 
as they cool and contract. Blue curves are isochrones at (bottom to top) 30, 160, and 300 
Myr; kinks in the older two isochrones arise from deuterium burning (objects burning D are 
substantially hotter than lower mass objects of the same age). Green curves are constant 
luminosity curves at (left to right) log L/L & = —5, —4.75, —4.5. For planet b solid contours 
denote fits to only the photometry while dashed curves are fits to photometry and H and 
K-band spectra. Crosses denote the individual model cases plotted in Figures 4-6. 
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Fig. 9. — Model spectra at fixed T e g = 900 K and varying gravities (labeled along right 
hand side), including several of the cases shown in Figure 2. Models are shown at a spectral 
resolution R = 1000. 
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Fig. 10. — Illustration of the effect of gravity and cloud properties on modeled emergent flux 
for T eS = 1000 K and log g = 5.0 (a) and 3.75 (b). Both plots (a) and (b) consist of four 
sub-panels. The right-most sub-panel depicts the T(P) profiles for three atmosphere models 
with the indicated T e g and log g. In both cases the profiles are for (left to right) for cloudless, 
/sed = 2, and 1 models. Thick lines denote the convective regions of the atmosphere models. 
The dotted line denotes chemical equilibrium between CO and CH 4 . The dashed lines are 
the condensation curves for Fe (right) and Mg 2 Si0 4 (left). The cloud base is expected at 
the point where the condensation curves cross the T(P) profiles. Remaining panels show the 
contribution function (see text) averaged over the J, H, K, L' and M' bandpasses (colored 
lines) for each of the three model cases. The shaded regions denote the extent of the cloud, 
extending from the point where the integrated optical depth from the top of the model is 
0.1 to the cloud base. Thick horizontal dashed line denotes cloud r = 2/3. 
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Fig. 11. — Definition of the transition from cloudy to cloudless surface boundary condition 
for the evolution. This represents a toy model of the L/T transition. In the hybrid toy model 



of Saumon & Marley (2008), the transition region is independent of gravity and the cloud 
clearing occurred between T e g = 1400 and 1200 K (lightly hashed area). To the right of 
the transition region shown, the surface boundary condition is based on cloudy atmosphere 
models; to the left, on cloudless atmospheres; and on a simple interpolation in the transition 
region. Here, we present an evolution calculation where the T c g range of the transition is 
made gravity dependent (densely hashed area). Representative cooling tracks are shown in 
black and labeled by the mass. Isochrones are the blue dotted lines. 



- 52 - 




Fig. 12. — Examples of cooling tracks for objects of 5Mj (red) and 20 Mj (blue) in a Mk vs. 
J — K (MKO system) color-magnitude diagram where the transition from cloudy (f se d = 1) 
to cloudless atmospheres is taken into account explicitly as in Saumon & Marley (2008). 
Dashed lines show the evolution when the transition occurs over a fixed range of T e g that 
is independent of gravity, solid lines show the evolution for the gravity- dependent transition 
(see Fig. 11). The planets in the HR 8799 planets are shown with green symbols while 



resolved field objects are shown in black (M dwarfs), red (L dwarfs) and blue (T dwarfs). 



The photometry is from 


Leggett et al. (2002), Knapp et al. (2004), Marocco et al. (2010) 


McCaughrean et al. 


2004 


, Burgasser et al. 


2006 


, and Liu & Leggett 


( 2005 ) . The parallaxes 


are from Perryman et al. 


1997), Dahn et al. 


(200 


2), Tinney et al. (2003), Vrba et al. (2004), 


Marocco et al. (201C 


), and various references in 1 


.eggett et al. (2002) 





